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1.0  SUMMARY 


Bodner's  constitutive  equation  [1-5]  is  one  of  those  used  by  the  U.S.  Air 
Force  to  model  viscoplastic  deformation.  This  particular  equation  is  useful 
because  of  its  relative  mathematical  simplicity,  partly  due  to  the  elimination 
of  the  need  for  a  yield  surface.  However,  Bodner's  formulation  is  not  without 
problems,  such  as  difficulties  in  fitting  the  material  constants  of  the  equa¬ 
tion  to  data  [6],  and  a  general  deficiency  of  the  model  to  properly  account  for 
hardening  phenomena  [7].  Deficiencies  in  the  hardening  behavior  of  Bodner's 
equation  likewise  imply  that  similar  deficiencies  probably  exist  for  the 
softening  behavior  caused  by  the  accumulation  of  material  damage  in  the  form  of 
voids  or  microcracks.  These  deficiencies  led  to  the  definition  of  the  present 
work,  the  objectives  of  which  were: 

(1)  to  improve  the  capabilities  of  Bodner's  viscoplastic  constitu¬ 
tive  equation  for  modeling  deformation-induced  hardening  phe¬ 
nomena  by  incorporating  microstructural  models  into  the  equa¬ 
tion,  and 

(2)  to  develop  a  continuum  formulation  that  simultaneously  repre¬ 
sents  the  mechanics  of  both  hardening  and  damage  behavior,  with 
the  ultimate  aim  of  developing  a  new  material  state  tensor. 

The  first  objective  has  been  accomplished.  The  second  objective  is  about 
50%  accomplished  at  this  writing,  and  work  required  for  completion  has  been 
defined.  These  results  are  briefly  described  in  the  following  paragraphs. 

Bodner's  constitutive  equation  in  uniaxial  form  is: 


where  e  is  the  strain  rate,  eQ  is  the  maximum  strain  rate,  C  and  n  are  material 
constants,  Z  is  the  hardening  parameter,  and  a  is  the  stress.  The  first 
objective  was  to  express  this  equation  in  a  form  that  would  readily  admit 


microstructural  models  from  the  metallurgical  literature.  This  was  accom¬ 
plished  using  some  recent  developments  in  mathematics  which  are  particularly 
suited  for  the  development  of  scaling  relations.  The  "new  mathematics"  are 
called  fractals,  and  are  geometric  forms  which  are  very  useful  for  concisely 
representing  the  complexity  of  nature.  The  reader  is  referred  to  [8]  and 
Section  2.2  of  this  report  for  descriptions  and  examples,  with  the  notation 
that  only  the  simplest  scaling  concepts  of  fractals  were  used  to  obtain  the 
following  result: 


e  exp 
o  h 


Physics  Data  „  ..  ^ 

Engrg.  Data 


(1.2) 


where  e$  is  the  macroscopic  strain  rate.  "Physics  Data"  are  known  quantities 
on  the  atomic  scale,  such  as  atomic  binding  energy.  "Engineering  Data"  are 
measured  quantities  on  the  macroscopic  scale,  such  as  stress.  M^,  Mo...  repre¬ 
sent  microstructural  phenomena  on  intermediate  scales,  described  by  models 
available  from  the  metallurgical  literature.  Examples  are  =  the  number  of 
atoms  per  unit  length  dislocation,  =  the  number  of  dislocations  per  disloca¬ 
tion  cell.  In  general,  any  number  of  M‘s  may  be  employed  for  arbitrarily 
selected  intermediate  scales  as  long  as  the  first  scale  relates  to  the  atomic 
behavior,  and  the  last  scale  relates  to  the  macroscopic  sample  behavior,  lhe 
formulation  is  quite  general.  Its  flexibility  permits  versatile  use  in  struc¬ 
tural  analysis  computer  codes  because  various  M's  may  be  inserted  in  or  deleted 
from  Equation  (1.2)  at  w.il  as  performance  variables  ( i . e . ,  temperature)  change 
during  a  simulation.  The  above  results  are  described  in  Section  3  of  this 
report. 

The  second  objective  was  to  develop  a  continuum  representation  that  unites 
deformation  (hardening)  and  damage  (softening)  mechanics,  with  the  ultimate  aim 
of  developing  a  new  material  state  tensor.  This  effort  began  with  an  analysis 
of  a  fractal  form  of  the  conservation  of  mass  equation,  but  the  first  attempt 
was  abandoned  when  it  was  found  that  hardening  and  damage  phenomena  remained 
separated  and  could  not  be  uniquely  combined  in  this  approach.  These  results 
are  described  in  more  detail  in  the  Appendix. 


The  second  attempt  to  accomplish  objective  (2)  was  successful.  This 
approach  combined  fractal  scalings  with  concepts  from  the  "nonlocal"  theories 
of  continuum  mechanics.  Nonlocal  theories  characterize  the  effects  of  material 
inhomogeneities  (i.e.,  hardening  and  damage  substructures)  on  macroscopic 
performance  by  including  a  fundamental  measure  of  length  in  the  formulations. 
This  length  parameter  characterizes  the  spacing  between  the  mi crostructural 
features  which  control  the  material  response,  or  the  ranges  of  the  fluctuations 
in  the  stress  and  strain  fields  within  the  material.  Internal  stress  field 
fluctuations  are  sometimes  called  "stress  polarization."  Insertion  of  fractal 
scaling  relations  into  Bodner's  Equation  resulted  in  the  development  of  a 
rather  complicated  expression  for  the  stress  polarization  (fluctuations),  ap, 
on  any  scale  m  within  the  deforming  material: 

apm  =  dEm/drm*  S’  °s  J  (U3) 

where  g(rm)  is  a  deformation-dependent  scale  function,  Em  is  the  effective 
modulus  at  scale  m<s ,  the  subscript  s  denotes  the  macroscopic  sample  scale,  and 
rm  is  a  function  of  the  volume  fractions  of  hardening  and  damage  substructures 
(•^hm.  Sim)  at  scale  m.  The  modulus  term  in  Equation  (1.3)  is  given  by  a  par¬ 
ticular  effective  medium  theory  [9,10],  Both  g(rm)  and  the  modulus  term  depend 
on  Sm»  Sm»  and  bm*  wbere  the  latter  is  a  numerical  scale  factor  (length 
parameter)  chosen  at  will.  Note  that  Sm»  Sm»  and  bm  are  a11  measurable,  so 
that  experimental  verification  is  possible.  The  function  rm  provides  the  basis 
for  defining  the  new  hardening/damage  state  tensor  in  terms  of  measures  of 
departure  from  an  identity  tensor,  which  in  turn  can  be  introduced  into  stan¬ 
dard  continuum  formulations.  Equation  (1.3)  thus  provides  the  foundation  for 
accomplishing  the  second  major  objective  stated  above. 

Because  of  the  mathematical  complexity  of  Equation  (1.3),  this  goal  could 
not  be  fully  attained  within  the  resources  of  this  project.  However,  prelimi¬ 
nary  results  concerning  the  numerical  verification  of  Equation  (1.3)  are  shown 
in  Section  4.2.  It  should  be  noted  that  accomplishment  of  the  second  objective 
will  also  provide  an  engineering  tool  which  apparently  does  not  presently 
exist:  a  nonlocal  theory  for  rate-dependent  deformation,  where  substructural 


"inhomogeneities"  evolve  rather  than  having  been  completely  specified  from  the 
initial  geometry  of  the  problem,  as  in  a  composite.  Engineering  applications 
of  such  a  development  include  prediction  of  the  deformation-dependent  spectrum 
of  scales  over  which  nonlocal  behavior  dominates.  This  in  turn  provides  guid¬ 
ance  concerning  the  proper  magnifications  needed  for  metal  1 ographi c  examina¬ 
tions  to  fully  reveal  the  substructural  features  that  control  deformation 
behavi or. 

The  balance  of  this  report  begins  with  the  Background,  which  includes 
brief  descriptions  of  Bodner's  Theory  and  fractals.  This  is  followed  by  a  sum¬ 
mary  of  the  work  performed  to  accomplish  the  first  objective  above:  the  incor¬ 
poration  of  microstructural  models  into  Bodner's  Equation.  Section  4  describes 
work  performed  toward  accomplishing  the  second  objective,  and  contains  the  com¬ 
plete  development  of  the  so-named  stress  polarization  equation.  Conclusions 
and  Recommendations  are  summarized  in  the  final  section. 


2.0  BACKGROUND 


2.1  BODNER'S  EQUATION  AND  THE  STATE  VARIABLE  APPROACH 

A  number  of  mathematical  approaches  a"e  available  for  the  constitutive 
modeling  of  rate-dependent  deformation  processes  [11],  each  having  unique 
advantages  and  disadvantages  concerning  physical  meaning,  developmental  costs, 
and  degree  of  success  in  application.  One  of  the  most  effective  approaches  is 
the  state  variable  method,  which  combines  macroscopic  stress-strai n  rate  rela¬ 
tions  with  judiciously  selected  microstructural  observations.  This  is  accom¬ 
plished  by  the  specification  of  both  a  kinetic  "equation  of  state"  which  con¬ 
tains  microstructural  "state  parameters",  and  supplementary  equations  describ¬ 
ing  tne  evolution  of  those  parameters  with  deformation.  An  example  of  such  a 
kinetic  equation  is  the  uniaxial  form  of  Bodner's  Law  shown  in  Equation  (1.1), 
which  describes  the  relation  between  plastic  strain  rate  ( e)  and  stress  (a). 
The  foundations  for  this  equation  are  summarized  in  the  following  two 
paragraphs. 

Briefly,  Bodner  decomposes  the  deformation  rate  tensor  into  elastic  and 
plastic  components  [1-5] 


(2.1) 


where  d^j  is  the  symmetric  part  of  the  material  velocity  gradient 
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The  flow  rule  employed  is  a  generalized  Prandtl -Reuss  law 


F(J2,Z,T)  sij 


(2.3) 
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where  S^j  is  the  deviatoric  stress  tensor,  J2  is  the  second  invariant  of  S^j,  T 
is  temperature,  and  Z  is  an  internal  state  parameter.  Squaring  Equation  (2.3), 


we  have 


F2  =  dP/J2 


where  0^  is  the  second  invariant  of  the  deformation  rate  tensor. 

A  key  feature  of  8odner's  theory  is  that  Of)  \ s  assumed  to  be  a  function  of 


02  =  f(J2) 


(2.5) 


so  that  the  above  F  can  be  defined  in  terms  of  stress.  The  specific  form 
chosen  for  Equation  (2.5)  was  motivated  by  dislocation  dynamics:  the  disloca¬ 
tion  velocity,  and  thus  the  strain  rate,  is  a  function  of  stress.  Equa¬ 
tion  (2.5)  can  be  stated  either  in  power  law  form  [12],  or  in  exponential  form 
[2,3]  as  follows: 


=  Dq  exp[-C(Z2/J2)n  j 


where  0Q,  C,  and  n  are  material  constants.  The  mathematical  simplicity  and 
clear  physical  meaning  of  Bodner's  model  permit  the  parameter  Z  to  be  inter¬ 
preted  as  a  "threshold  stress"  or  as  a  hardening  parameter  in  extensions  of  the 
theory  to  strain  rate  and  strain  rate  history  dependencies  [12,13].  In  recent 
papers  [3,4],  Bodner  defined  Z  as  a  hardening  parameter  describing  "the  micro- 
structural  state  of  the  material  at  time  t"  or  the  "stored  energy  of  cold 
work",  and  extended  the  formulation  to  anisotropic  conditions  by  considering  Z 
to  be  a  tensor. 

However,  some  problems  still  remained.  It  was  observed  that  correlations 
for  Z  were  apparently  not  transferable  from  one  data  set  to  another,  even  for 
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the  same  material.  Additionally,  some  hardening  phenomena  were  not  properly 
modeled  by  Equation  (2.6),  or  Equation  (1.1).  This  may  have  been  caused  by 
Bodner's  use  of  the  plastic  work  as  a  measure  of  the  viscoplastic  process 
[7,14],  However,  it  also  seemed  entirely  possible  that  the  description  of  Z 
was  simply  incomplete,  as  explained  in  the  following  paragraphs. 

Recent  conferences  such  as  [15]  have  recognized  the  importance  of  includ¬ 
ing  mathematical  representati ons  of  mi crostructural  phenomena  in  macroscopic 
constitutive  equations.  In  particular,  see  the  articles  by  Drucker,  Hart,  and 
Aifantis  in  [15].  However,  an  approach  used  in  many  current  constitutive 
models  is  to  minimize  the  number  of  internal  state  parameters.  The  reasons  for 
this  approach  are  that  (1)  larger  numbers  of  parameters  can  lead  to  non-unique 
parameter  values  in  benchmarking  exercises,  and  (2)  physical  justifications  are 
increasingly  harder  to  find  as  the  number  of  mathematical  "free  parameters" 
increases.  Although  Bodner  has  recently  developed  a  tensor  form  of  the  harden¬ 
ing  parameter  Z,  each  element  of  the  tensor  can  be  interpreted  as  a  single 
macroscopic  state  parameter  expressed  in  terms  of  five  empirically  determined 
macroscopic  material  constants. 

Consequently,  some  investigators  have  viewed  the  state  variable  method  as 
..cving  "reached  a  point  of  diminishing  returns"  [16]  because  of  the  increasing 
mathematical  complexity  required  for  improved  predictive  accuracy,  and  the  lack 
of  agreement  concerning  which  state  parameters  are  physically  important  and 
computationally  useful.  These  difficulties  arise  because  the  various  micro- 
structural  phenomena  characteri zed  by  state  parameters  may  each  dominate  the 
flow  process  at  different  stages  of  deformation  or  in  different  temperature 
regimes.  This  in  turn  may  require  that  a  new  constitutive  equation  be 
developed  for  each  regime  of  performance  for  each  material,  even  when  micro- 
structural  correlations  are  available  over  these  regimes.  The  resulting  mathe¬ 
matical  complexity  can  indeed  be  discouraging. 

Section  3  of  this  report  describes  an  alternative  approach  to  the  formula¬ 
tion  of  such  constitutive  equations,  with  the  intent  of  providing  a  method  that 
will  reduce  the  recurring  developmental  requirements  described  above,  in  addi¬ 
tion  to  possibly  improving  the  understanding  of  how  mi crostructural  evolution 
determines  macrostructura 1  response.  These  developments  rely  on  a  recently 


developed  form  of  geometry,  known  as  fractals.  A  brief  summary  of  the  attri¬ 
butes  of  fractals  which  are  of  use  in  this  work  is  provided  in  the  following 
paragraphs . 

2.2  FRACTALS 

A  fractal  is  a  geometric  form  which  possesses  two  important  and  useful 
attributes:  noninteger  dimensional ity  and  scale  invariance.  These  attributes 
may  be  best  described  by  the  examples  shown  in  Figures  2.1  and  2.2.  Figure  2.1 
shows  a  fractal  known  as  a  Koch  Curve  as  it  develops  during  the  iterative  con¬ 
struction  process.  Beginning  with  a  line  of  unit  length  (Figure  2.1a),  the 
line  is  divided  into  b  segments  of  equal  length,  and  then  the  original  line  is 
replaced  by  N  new  segments  of  length  1/b,  as  shown  in  Figure  2.1b.  The  next 
iteration  is  shown  in  Figure  2.1c,  where  each  segment  is  again  subdivided  and 
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FIGURE  2.1.  A  Koch  Curve  (0=1.26) 


replaced  by  N  new  segments.  As  this  process  continues  indefinitely,  the  frac¬ 
tal  develops  scale  invariance.  That  is,  the  fractal  has  the  same  appearance 
under  any  magnification. 

The  fractal  dimension  is  defined  rigorously  in  [8]  using  topological 
arguments.  For  the  present  purposes,  computation  of  the  fractal  dimension  can 
be  briefly  described  as  follows: 


b°  =  N 


This  means  that  the  basic  scaling  unit  raised  to  the  power  of  the  dimension  is 
equal  to  the  number  of  new  units  (segments).  In  Figure  2.1,  b=3  and  N=4,  so 
that  the  fractal  dimension  of  that  particular  curve  is  D  =  ln(4)/ln(3)  =  1.26. 
The  result  is  that  the  Koch  Curve  is  neither  a  line  nor  a  plane,  but  is  some¬ 
where  in  between.  In  essence,  this  fractal  dimension  describes  how  much  space 
the  curve  fills  in  following  its  meandering  course. 

A  second  example,  known  as  a  Sierpinski  Carpet,  is  shown  in  Figure  2.2. 

In  this  case,  a  square  is  divided  into  b=3  segments  on  each  side,  resulting  in 
nine  subsquares.  Then  1=1  subsquare  is  removed.  The  removed  subsquare  is 
called  a  "tremma",  and  it  creates  empty  space  in  the  fractal,  which  is  called 
"lacunarity"  [8],  a  term  apparently  originating  from  the  description  of  cloud 
structures  in  meteorology.  Figure  2.2b  shows  the  Sierpinski  Carpet  at  the 
second  stage  of  construction.  When  this  process  is  repeated  indefinitely,  the 
fractal  again  becomes  scale  invariant.  In  this  case,  the  fractal  dimension 

D  ?  9 

must  be  clarified:  b  =  the  number  of  remaining  new  units,  which  is  b  -  1  , 
so  the  fractal  dimension  of  this  particular  Sierpinski  Carpet  is 


0  =  1 n (b2-l 2  )/l n( b)  =  1.89 


Three  items  are  worth  noting.  First,  fractals  can  be  easily  extended  to 
higher  dimensions.  For  example,  the  Sierpinski  Carpet  embedded  in  Euclidean 
3-space  becomes  a  Sierpinski  Sponge  [8]  with  dimension  between  2  and  3. 
Second,  for  the  purposes  of  this  work,  it  is  convenient  to  recognize  that 


fractals  can  be  constructed  by  either  adding  or  subtracting  substructures, 
which  results  in  the  fractal  dimension  being  either  greater  than  or  less  than 
the  dimension  of  the  Euclidean  space  in  which  it  is  embedded.  This  feature 
will  be  of  use  in  later  discussions  concerning  the  differences  between 
deformation  and  damage.  The  third  item  is  a  matter  of  nomenclature.  When  a 
fractal  is  iterated  a  finite  number  of  times,  it  is  called  a  fractal  lattice. 
On  such  a  lattice,  scale  invariance  is  valid  only  between  the  two  extreme 
characteristic  lengths,  and  the  fractal  dimension  is  approximate  [8]. 

The  use  of  fractal  concepts  to  incorporate  mi crostructural  models  into 
Bodner's  constitutive  equation  is  described  in  the  next  section. 
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3.0  A  MICROSTRUCTURAL  FORM  OF  BODNER'S  EQUATION 


The  objective  of  this  section  is  to  develop  and  verify  a  new  form  of 
Bodner's  Equation  with  improved  capabilities  for  modeling  the  hardening 
behavior  of  viscoplastical ly  deforming  materials.  It  was  believed  that  such 
improvements  could  be  obtained  by  incorporating  mi crostructural  information 
into  Bodner's  constitutive  equation.  In  this  manner,  Bodner's  continuum  repre¬ 
sentation  could  benefit  from  the  numerous  results  of  metallurgical  research. 
Because  of  the  difficulties  associated  with  the  state  variable  method  (see 
Section  2,1),  an  alternative  approach  was  taken  to  accomplish  this  task.  The 
basic  concept  was  as  follows.  It  has  been  observed  that  new  substructures  of 
(apparently)  increasing  complexity  are  created  as  materials  deform  visco- 
plastically.  This  is  somewhat  analogous  to  the  construction  process  for  a 
fractal,  as  outlined  in  Section  2.2.  It  was  believed  that  the  scaling  proper¬ 
ties  (self-similarity)  of  fractals  could  be  advantageously  employed  to  relate 
the  observed  mi crostructural  forms  to  macroscopic  performance  variables.  Per- 
suit  of  this  approach  led  to  the  replacement  of  Bodner's  empirical  Z  and  n  by 
known  metallurgical  formulas  describing  mi crostructural  evolution.  This  new 
constitutive  form  was  expected  to  exhibit  improved  hardening  capabilities,  thus 
accomplishing  the  above  objective. 

The  outline  of  this  section  is  as  follows.  The  hardness/stress  ratio  of 
Bodner's  Equation  was  interpreted  as  a  ratio  of  potential  gradients  per  unit 
volume.  Fractal  scaling  concepts  were  then  employed  to  render  the  stress  equal 
on  all  scales,  thus  establishing  local  behavior  with  respect  to  scale.  This  in 
turn  permitted  an  expansion  of  the  above  potential  ratio,  resulting  in  Bodner's 
Equation  expressed  as  a  product  of  microstructural  flow  resistance  contribu¬ 
tions.  An  approximation  of  this  expanded  form  in  terms  of  fractal  structures 
then  produced  a  constitutive  equation  explicit  in  scale-dependent  microstruc¬ 
tural  quantities  such  as  numbers  or  densities  of  atomic  bonds,  dislocations, 
cells,  and  grains.  The  above  arguments  also  revealed  that  the  power  law 
exponent,  usually  found  empirically,  could  be  defined  from  products  of  known 
microstructural  models. 
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An  example  application  of  this  new  constitutive  form  to  the  creep  of  a 
20%Cr-35%Ni  steel  alloy  provided  limited  proof  that  the  methodology  was 
valid. 

3.1  DEVELOPMENT 

The  constitutive  developments  begin  with  a  particular  form  of  Bodner's 
Equation  that  separates  out  the  hardness/stress  ratio.  For  uniaxial  loading, 
this  form  is: 


-C  «n(es/e0)  =  (Z/^)' 


(3.1) 


where  =  macroscopic  plastic  strain  rate,  =  limiting  strain  rate  in  shear, 
C  =  a  constant,  as  =  macroscopic  stress,  Z  is  a  hardness  parameter  representing 
the  "stored  energy  of  cold  work"  [4,5],  and  n  is  an  empirically  determined 
exponent.  Note  that  the  left  side  of  Equation  (3.1)  is  a  representation  of  a 
flow  rate  ratio  and  the  right  side  is  a  representation  of  a  potential  ratio. 
Both  Z  and  a$  have  units  of  stress,  which  is  energy  E  (or  potential)  per  unit 
volume  V.  Flow  rates  are  frequently  expressed  as  functions  of  potential 
gradients  [17,18],  and  it  is  postulated  that  the  right  side  of  Equation  (3.1) 


can  be  written  as: 


(kf  *  f 


where  cra  is  a  potential  per  unit  volume  on  the  atomic  scale  because  all 
deformation  ultimately  involves  atomic  bond  stretching  or  breaking. 

For  each  scale  i,  the  stress  is 


TO 


L.j  is  a  characteristic  length  and  is  the  fractal  dimensionality.  The  frac¬ 
tal  volumes  serve  as  scaling  factors  for  the  stresses.  Proper  (but  unknown) 
choices  of  will  render  the  stress  invariant  with  respect  to  scale.  That  is, 

any  arbitrarily  chosen  subsample  of  size  between  La  and  Ls  will  have  the  same 
stress,  as  in  a  true  continuum.  Note  that  although  a.j/Cj  =  1  for  any  i,j,  it 


is  not  necessarily  true  that  Sk/o^  =  1. 


The  above  condition  essentially  establishes  local  behavior  with  respect  to 
scale,  similar  to  local  behavior  with  respect  to  space  [19-22].  Mathemati¬ 
cally,  scale  invariant  local  behavior  can  be  expressed  as  (see  Discussion) 


a  =f(a(a(a(a)))) 
a  dv  cv  g  s' ' 


(3.5) 


where  d,c,g  are  the  dislocation,  cell  and  grain  scales,  respectively.  Other 
scales  are  also  possible.  Then  applying  the  chain  rule  to  Equation  (3.2) 
yields  the  product  form: 


5a  da  5a  5a 

-C  *n(e s/e  )  =  ^  ^  =  n  (Wl} 

d  c  g  s  i 


(3.6) 


where  i  =  a,d,c,g  and  n  is  the  product  operator.  Nonlocal  behavior  would  be 
manifested  as  additive  terms  in  Equation  (3.6)  via  the  chain  rule. 

The  last  term  of  Equation  (3.6)  can  be  expanded  using  Equation  (3.3): 


Each  of  the  in  Equation  (3.7)  contains  a  fractal  dimensional ity  D^.  It  is 
natural  to  assume  that  D  changes  with  both  scale  and  deformation.  However,  the 
D's  are  assumed  unknown  in  the  present  development,  and  will  be  eliminated  with 
the  approximation  shown  by  the  dashed  lines  in  Figure  3.1,  as  follows.  First, 
substitution  of  Equations  (3.7)  into  (3.6)  yields,  after  cancellations, 


where  i  =  a,d,c,g,  and  the  subscript  j  is  defined  in  Figure  3.1.  Note  that 
Equation  (3.8)  contains  products  of  lengths  of  powers  (local),  rather  than  a 
sum  of  lengths  to  powers  (nonlocal)  as  suggested  in  [23].  Second,  the  charac¬ 
teristic  lengths  may  be  expressed  in  a  standard  form  by  scaling  upward  from 
the  characteristic  length  of  the  atomic  scale,  La  =  a0,  an  appropriate  lattice 
constant.  Using  the  scaling  parameter  (b)  of  Equation  (2.7),  we  have 
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FIGURE  3.1.  Conceptual  Diagram  for  Dimensional ity  as  a  Function  of  Scale 
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(3.9) 
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where  m^  is  the  number  of  fractal  construction  iterations  required  to  scale 
from  the  atomic  to  the  ith  scale.  Consequently,  the  length  ratios  in  Equa¬ 
tion  (3.8)  can  be  written 


D  .  Am .  D  . 

(WLi>  J  -  b  J  J 


(3.10) 


where  Amj  is  the  number  of  iterations  between  scales  i  and  i+1.  Lastly,  sub¬ 
stituting  Equations  (2.7)  and  (3.10)  into  (3.8)  eliminates  the  D ' s  if  b  is  held 
constant  over  all  scales,  yielding 

.  .  /£a  -  V  \  Am, 

-e  *i(,  /c  )  ■(  - - 5^)  n  (N.  J)  (3.11) 

VES  -  v's/  j  J 

Am . 

where  N-  J  is  the  total  number  of  new  fractal  structure  subunits  introduced  in 

J 

iterating  from  scale  i  to  i+1.  These  N-terms  represent  microstructural  fea¬ 
tures  such  as  the  number  of  atoms  per  dislocation  or  the  number  of  dislocations 
per  dislocation  cel  1 . 

Additional  approximations  reduce  this  constitutive  equation  to  a  simpler 
form.  The  terms  Vi  represent  the  rates  of  volumetric  change,  and  can  be  set  to 
zero  by  virtue  of  the  first  strain  invariant  (^=0),  a  frequent  assumption 

employed  in  non-elastic  analyses.  Expressing  Es  as  the  plastic  work  yields 

•  • 

Es  =  ^ s  es,  where  =  <j>cxs  and  <t>  is  the  Schmid  factor.  On  the  atomic  scale, 

• 

Ea  =  Fa  va,  and  the  separation  force  can  be  approximated  as  Fa  =  AEa/Axa  = 
E0/a0,  where  EQ  is  the  atomic  binding  energy.  The  velocity  va  can  be  assumed 
equal  to  the  average  velocity  of  mobile  dislocations  vm.  Substitution  of  these 
approximations  into  Equation  (3.11)  yields,  after  rearranging: 
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e  exp 
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n  (N-  J) 
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3.5 


(3.12) 


The  exponent  in  Equation  (3.12)  has  essentially  been  separated  into  the 
form:  (Physics  Data)/Engineering  Data)  x  (Metallurgical  Models),  where 
"Physics  Data"  represents  known  quantities  on  the  atomic  scale,  "Engineering 
Data"  represents  known  measurables  on  the  macroscopic  sample  scale,  and  "Metal¬ 
lurgical  Models"  represents  microstructural  phenomena  on  intermediate  scales. 

If  the  mi crostructural  phenomena  controlling  the  deformation  rate  are  known, 
models  from  the  metallurgical  literature  may  be  substituted  for  the  N-terms. 
Different  substitutions  can  be  made  as  the  controlling  microstructure  evolves 
with  deformation  or  temperature.  Equation  (3.12)  thus  provides  a  flexible 
methodology  for  incorporating  microstructural  information  into  macroscopic 
constitutive  equations. 

3.2  APPLICATION 

Limited  verification  of  Equation  (3.12)  was  obtained  by  simulating  the 
hardening  behavior  for  a  20%Cr-35%Ni  steel  undergoing  a  normal  creep  transient 
and  entering  the  steady  state  creep  regime  at  700°C  [24,25,26],  as  follows. 
First,  the  dislocation  velocity  vm  in  Equation  (3.12)  was  related  to  the  macro- 
scopic  strain  rate  by  the  Orowan  equation,  =  Bpmvm,  where  B  is  the  Burger's 
vector  and  pm  is  the  density  of  mobile  dislocations  (cm-4).  This  cancelled 
from  the  right  side  of  Equation  (3.12).  Only  atomic  and  dislocation  N-terms 
were  employed  in  this  first  approximation,  and  both  were  assumed  to  represent 
stored  energies,  consistent  with  the  definition  of  Z  [5],  N^  represented  the 
number  of  atomic  bonds  left  in  a  strained  condition  after  a  dislocation  loop 
expands  to  a  typical  radius  proportional  to  the  total  dislocation  density,  r  = 
PT'1/2  [25]: 

N^1  =  2n/(aQPT1/2)  (3.13) 

The  second  microstructural  term  represented  the  energy  storage  at  the  disloca¬ 
tion  scale,  and  was  defined  as  the  number  of  immobile  dislocations  per  square 
cm  of  sample  area: 


3.6 


3. 


Substituting  the  above  microstructural  models  into  Equation  (3.12)  gave 


V  eo  exp 


2uCEopI 

2  1/2 
a„B<0a  p  pT1,£ 

o  s  m  T  > 


(3.15) 


The  constitutive  equation  was  then  expressed  in  terms  of  only  the  total 
dislocation  density  Py,  a  more  desirable  state  parameter  [11],  by  employing  the 
dislocation  distribution  models  in  [27]: 


Pn/PT 


=  exp(-kQ)  y  * 


k-r) ! 


(3.16) 


PI  =  PT  “  pm 


(3.17) 


where  9  =  GBpy1/,2/Ts,  G  is  the  shear  modulus,  and  k  is  an  integer  describing 
the  "sharpness"  of  the  distribution  [27].  Substituting  Equations  (3.16)  and 
(3.17)  into  (3.15)  gave 


exp  - 


aoBpas PT 


o  1-f 
~T/2  f 


(3.18) 


where  f  is  defined  by  the  right  hand  side  of  Equation  (3.16).  Estimated  values 
for  the  above  material  parameters  are  [25,26,28]:  eQ  =  3.13  x  10-5/sec,  E0  = 
4.29  eV/atom  at  0°K,  a0  =  2.86  x  10"B  cm,  B  =  2.5  x  10  B  cm,  G  =  7.3  x 

c  o 

10  kg/cm  ,  i>  =  1/3.1.  The  constant  C  was  taken  as  unity,  and  k  in  Equa¬ 
tion  (3.16)  was  set  equal  to  4  based  on  observation  of  the  dislocation  distri¬ 
butions  in  [26]  compared  to  those  in  [27]. 
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;uation  (3.18)  was  tested  in  the  following  manner.  Th< 
ir  Fe-20%/Cr-35foNi  at  700°C  and  1300  kg/cm^  was  determ' 
g.  2  of  [25],  and  was  used  as  input  to  solve  iterativf 
g  Py  versus  time.  The  results  are  shown  in  Figure  3.; 
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another  prediction  from  [25],  and  data  at  an  applied  stress  of  1100  kg/cm2 
[26],  Unfortunately,  data  at  1300  kg/cm2  were  not  available  for  comparison. 
However,  Figure  3.2  shows  reasonably  good  quantitative  agreement  with  the  pre- 

o 

diction  of  [25]  at  1300  kg/cnr,  which  required  three  adjustable  parameters  and 
reproduced  the  creep  history  rather  well.  Although  the  qualitative  agreement 
with  the  data  at  the  lower  stress  is  acceptable,  improvements  are  possible  by 
decreasing  k  and  adjusting  Py  as  the  time  approaches  zero.  The  "softer" 
distribution  at  lower  strains  (shorter  times)  is  also  documented  in  [26]. 

The  above  results  are  subject  to  errors  originating  from  the  graphical 
interpretation  of  cs.  approximati on  in  Figure  3.1,  and  the  approximation 
that  mi crostructural  phenomena  on  scales  ’arger  than  dislocations  were  negli¬ 
gible.  However,  the  reasonably  good  results  can  be  taken  as  an  indication  that 
the  above  developments  are  valid  and  useful  constitutive  forms. 

3.3  DISCUSSION 

The  state  variable  approach  to  developing  viscoplastic  constitutive  equa¬ 
tions  is  a  valuable  method  because  it  has  the  capability  of  combining  models 
from  both  continuum  mechanics  and  metallurgy,  resulting  in  a  pragmatic  use  of 
available  information.  However,  the  resulting  complexity  of  the  mathematical 
expressions  is  sometimes  a  disadvantage,  and  the  physical  meaning  of  the  cor¬ 
relations  is  not  always  clear.  This  section  has  attempted  to  address  these 
problems  by  restating  Bodner's  macroscopic  constitutive  equation  in  a  form  that 
readily  admits  mi crostructural  concepts. 

A  number  of  qualifying  statements  are  in  order.  First,  the  power  law  form 
chosen  as  the  starting  point  of  this  development  [Equation  (1.1)]  limits  the 
present  results  to  uniaxial  loading  conditions.  Multiaxial  and  anisotropic 
forms  of  Bodner's  law  are  available  in  [4],  The  hardness/stress  ratio  of  Equa¬ 
tion  (3.1)  results  in  the  N-terms  being  defined  as  flow  resistance  or  energy 
storage  contributions,  whereas  other  ratios  may  require  different  definitions. 

Second,  there  is  currently  no  proof  that  Equation  (3.2)  is  true,  although 
the  validity  of  the  resulting  constitutive  development  is  supported  by  the 
example  application  described  above.  Further  testing  is  needed  for  a  variety 
of  materials  and  conditions.  Additional  work  is  also  needed  to  confirm  the 


following  observation:  the  expanded  form  in  Equation  (3.6)  implies  that  the 
empirical  power  law  exponent  (n)  effectively  originates  from  the  product  of 
microstructural  flow  resistance  contributions  on  various  scales  within  the 
material.  A  physical  interpretation  is  that  the  N-terms  (microstructural 
models)  are  commonly  expressed  as  functions  of  os  to  rational  powers  p  in  the 
metallurgical  literature.  After  taking  the  derivatives  in  Equation  (3.6),  n 
effectively  becomes  a  sum  of  p's.  Because  various  mi crostructural  phenomena 
dominate  at  various  stages  of  deformation,  this  could  possibly  explain  why  n 
frequently  assumes  non-interger  values  that  change  with  deformation.  Nonlocal 
behavior  could  result  in  a  more  complicated  n  =  f(p). 

Third,  the  time  derivatives  in  Equation  (3.2)  may  raise  questions  related 
to  Prager's  consistency  condition  [19].  This  form  was  postulated  in  order  to 
eliminate  the  plastic  work  integral,  which  is  an  "undesirable"  state  variable 
[7,11].  Although  other  gradients  may  be  possible,  the  present  approach  facili¬ 
tates  simpl if ication  of  the  right  hand  side  of  Equation  (3.11)  by  metallurgical 
rate  laws  and  the  conservation  of  volume. 

Fourth,  many  fractals  lack  translational  invariance  (i.e.,  p.  426  of  [8]), 
which  conflicts  with  the  fundamental  postulate  of  frame  indifference  in  con¬ 
tinuum  mechanics.  However,  the  fractal  scalings  render  the  material  behavior 
local  with  respect  to  scale,  thus  imposing  a  condition  similar  to  a  "true  con¬ 
tinuum"  (i.e.,  a  fluid)  which  possesses  such  invariance.  Additionally,  many 
useful  constitutive  theories  violate  this  postulate  [19],  so  the  question  may 
be  irrevelant  in  practice. 

This  problem  can  also  be  addressed  by  an  alternative  approach  which  has 
been  successful  in  solid  state  physics  [i.e.,  29],  That  approach  is  to  con¬ 
sider  the  viscoplastic  deformation  process  to  occur  on  fractal  lattices,  or  in 
"fractal  space,"  as  follows.  The  potentials  (E)  in  Equation  (3.3)  contain  only 
the  fundamental  quantities  of  force  and  length,  which  are  invariant  with 
respect  to  transf ormati ons  from  Euclidean  to  fractal  space.  However,  the 
volume  terms  contain  the  dimens ional i ty  (D),  which  can  only  assume  integer 
values  in  Euclidean  space.  Recognition  of  D  as  another  fundamental  quantity 
that  can  assume  non-integer  values  requires  the  proper  space  for  representa¬ 
tion.  The  above  constitutive  form  can  be  developed  by  employing  fractal  space. 


and  assuming  that  Z/cr  appears  as  a  potential  gradient  there,  as  in  the  central 
terms  of  Equation  (3.6).  Interpretation  of  these  gradients  (da^/^^j)  as  mea¬ 
sures  of  energy  transfer  from  one  scale  to  the  next  naturally  leads  one  to 
expect  local  behavior  with  respect  to  scale,  because  relations  between  remote 
scales  would  not  be  expected  to  be  significant  on  physical  grounds.  For 
example,  atomic  bond  breaking  would  have  a  greater  effect  on  dislocation 
dynamics  than  on  grain  size.  In  any  case,  the  number  of  terms  in  the  chain 
rule  expansion  can  easily  be  chosen  to  enforce  local  behavior.  The  result  is 
that  the  chain  rule,  which  is  still  valid  in  fractal  space  [29],  again  produces 
the  product  form  of  Equation  (3.6).  The  key  assumption  concerns  the  existence 
of  the  potential  gradient  in  fractal  space,  and  the  key  result  is  the  "con- 
tinuization"  of  the  stress  to  produce  local  behavior  as  viewed  from  fractal 
space. 

In  conclusion,  the  above  constitutive  form  does  not  provide  any  new 
information  about  mi crostructural  models  per  se.  However,  it  does  appear  to 
provide  a  useful  methodology  for  i ncorporati ng  those  models  into  macroscopic 
constitutive  equations,  and  a  physical  description  for  the  origin  of  the  power 
law  exponent.  The  results  of  these  developments  are  expected  to  be  manifested 
as  improvements  in  the  ability  of  Bodner's  constitutive  equation  to  model  the 
hardening  behavior  of  engineering  metals  and  alloys. 


4.0  PROGRESS  TOWARD  A  COMBINED  HARDENING/DAMAGE  STATE  TENSOR 


The  developments  in  Section  3  employed  an  approximation  of  the  scale- 
dependent  fractal  dimensionality  (D)  (Figure  3.1)  in  order  to  express  Bodner's 
Equation  in  terms  of  mi crostructural  models.  Although  the  above  results  are 
useful,  as  outlined  in  the  Summary,  it  was  thought  that  a  constitutive  equation 
addressing  both  hardening  and  damage  could  be  obtained  without  having  to 
incorporate  specific  metallurgical  models.  In  essence,  the  mi crostructural 
phenomena  would  be  implicit  in  the  constitutive  formulation,  which  in  turn 
would  be  a  function  of  only  well-known  data  on  the  atomic  scale  (i.e.,  atomic 
binding  energy)  and  the  macroscopic  state  variables  (i.e.,  strain  and  stress). 
The  initial  point  of  departure  for  this  attempt  was  Bodner's  Equation  in  the 
form  (see  Equation  (3.8)). 
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With  characteristic  lengths  and  energies  comparatively  well  known,  the  problem 
was  to  find  solutions  for  the  dimensionalities  in  Equation  (4.1)  in  terms  of 
only  the  macrostructural  state  variables  (es,  e$,  a5).  If  increases  or 
decreases  in  the  D‘s  were  assumed  to  represent  hardening  or  damage,  respec¬ 
tively,  these  solutions  could  then  be  used  to  define  a  combined  hardening/ 
damage  state  tensor.  An  unsuccessful  attempt  to  find  these  solutions  for  the 
uniaxial  isotropic  case  by  analysis  of  the  conservation  of  mass  equation  is 
described  briefly  in  the  Appendix. 

A  second  attempt  was  successful.  Concepts  from  the  noncl assical  ,  or  "non¬ 
local,"  continuum  theories  [20-22,  30-35]  were  found  to  provide  a  framework 
useful  to  accomplish  the  second  objective  in  the  Summary.  The  central  idea  was 
that  the  hardening  and  damage  substructures  act  as  inhomogeneities  that  cause 
fluctuations  in  the  internal  stress  fields,  or  the  "stress  pol ari zat i on" 
depicted  in  Figure  4.1.  The  problem  was  thus  redefined  to  be  the  development 


FIGURE  4.1.  Stress  Polarization  as  a  Function  of  Scale 


of  a  scale-dependent  viscoplastic  stress  polarization  equation  and  its  solution 
as  a  function  of  suitable  hardening  and  damage  state  parameters. 

The  general  approach  was  as  follows.  As  in  Section  3,  the  hardness/stress 
ratio  was  represented  by  a  ratio  of  potential  gradients,  except  that  nonlocal 
behavior  was  permitted  by  letting  the  numerator  represent  the  scale  at  which 
the  microstructure  controls  the  strain  rate.  This  was  used  to  develop  a  rela¬ 
tion  between  the  energy  storage  rates  for  the  controlling  and  macroscopic 
scales.  Fractal  scaling  concepts  were  then  employed  to  define  the  scale  factor 
for  the  above  relation  in  terms  of  hardening  and  damage.  Definition  of  a 
strain  rate  scaling  equation  then  resulted  in  scale  invariance  of  the  energy 
storage  rate,  with  respect  to  the  strain  on  each  scale.  Stored  energies  were 
defined  in  terms  of  effective  moduli  on  each  scale,  and  the  stress  polarization 
equation  was  produced  by  rearrangement  of  the  above  scale  invariant  storage 
rate  equation.  Embedded  within  the  polarization  equation  is  a  function 
describing  the  combined  hardening/damage  state  in  terms  of  departures  from 
unity,  thus  providing  the  basis  for  the  desired  new  state  tensor.  The  above 
developments  are  described  in  more  detail  in  the  next  subsection,  followed  by 
preliminary  verification  results. 


4.2 


4.1  A  STRESS  POLARIZATION  EQUATION  FOR  VISCOPLASTICITY 


A  key  requirement  for  expression  of  the  stress  polarization  equation  is 
the  following  scale  invariant  energy  storage  rate  equation 
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Although  this  equation  could  probably  have  been  simply  assumed  by  virtue  of  the 
“fractal  space"  arguments  in  Section  3.3,  the  following  derivation  is  presented 
to  place  the  polarization  equation  on  a  firmer  foundation. 


4.1.1  An  Energy  Storage  Rate  Scaling  Equation 


Arguments  begin  from  the  same  starting  point  as  did  the  microstructura 1 
developments  in  Section  3  of  this  report.  Bodner's  Equation  is  stated  in  the 
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where  the  subscript  s  denotes  the  sample  scale,  Z  is  the  hardening  parameter 
which  represents  the  storage  of  cold  work,  and  n  is  the  empirical  exponent.  It 
is  again  assumed  that  the  hardness/stress  ratio  on  the  right  hand  side  of  the 
above  equation  can  be  expressed  as  a  ratio  of  potential  gradients: 
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However,  in  this  case  it  is  assumed  that  the  potential  occurs  on  the  scale  m 
at  which  the  microstructure  is  controlling  the  deformation  rate  c  at  a  par¬ 
ticular  macroscopic  state  of  the  material  (es,  as).  In  essence,  m  represents 
the  scale  at  which  a  nonlocal  theory  is  required  to  describe  the  internal 
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stress  field  fluctuations  caused  by  hardening  or  damage  inhomogeneties.  Above 
this  scale,  local  theories  would  suffice.  It  is  assumed  that  the  scale  m  can 
change  as  the  material  state  changes,  i.e.,  as  deformation  advances. 


As  before,  each  potential  is  expressed  as  the  ratio  of  a  characteristic 
energy  over  a  characteristic  volume 
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where  the  volume  on  scale  i  is  a  fractal  volume  denoted  by 
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and  Li  is  the  characteristic  length  of  scale  i.  Substituting  Equation  (4.5) 
into  Equation  (4.4)  and  performing  the  differentiations.  Equation  (4.3)  can  be 
written  as 
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where  the  (*)  denotes  time  differentiation.  Note  that  because  the  chain  rule 
was  not  employed  and  differentiation  was  performed  independently  for  the 
numerator  and  denominator,  it  is  not  necessary  to  assume  local  behavior  with 
respect  to  scale. 

The  term  Em  in  the  numerator  represents  the  rate  of  energy  storage  by 

microstructural  features  on  scale  m,  a  generalization  of  the  results  in  [34], 

•  • 

and  is  hereafter  denoted  Um  for  clarity.  The  Es  term  in  the  denominator 
represents  the  rate  at  which  plastic  work  is  performed  by  the  effective  stress 
'F 0  ,  where  Y  denotes  the  Schmid  factor.  The  denominator  is  expressed  hereafter 


As  before,  the  first  strain  invariant  is  assumed  to  be  zero  (Tr  e  =  0)  on 
all  scales,  which  results  in  volume  conservation,  so  that  the  V  terms  in 
Equation  (4.7)  vanish.  When  the  above  substitutions  and  Equation  (4.6)  are 
inserted  into  Equation  (4.7),  the  result  after  some  minor  manipulations  is 
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This  equation  is  interpreted  to  mean  that  the  rate  of  plastic  work  on  the 
macroscopic  scale  (s)  times  a  rate  factor  equals  the  rate  of  energy  storage  on 
another  scale  (m<s)  times  a  scale  factor.  In  essence,  the  same  result  found  in 
[34]  has  been  expanded  to  include  both  rate  and  scale  dependence.  This  is 
valid  for  any  scale,  so  setting  m=s  in  Equation  (4.8)  gives 
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which  means  that  the  rate  of  macroscopic  plastic  work  times  the  rate  factor 
equals  the  rate  of  macroscopic  energy  storage.  This  result  seems  reasonable, 
so  that  combining  the  last  two  equations  produces 
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which  relates  the  macroscopic  energy  storage  rate  to  that  on  any  other  scale. 
This  equation  indicates  that  energy  storage  rate  (flow  resistance)  contribu¬ 
tions  from  different  scales  are  indeed  different,  and  require  quantification  by 
the  N-terms  in  Section  3  when  incorporated  into  Bodner's  Equation.  When  the  e 
of  Equations  (4.8)  and  (4.9)  is  retained,  a  more  useful  form  of  Equation  (4.10) 
is  produced 


(4.11) 
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4.1.2  Definition  of  the  Scale  Factors 

The  next  step  in  the  development  of  the  scale  and  rate  dependent  stress 
polarization  equation  concerns  the  definition  of  the  scaling  factor  in  paren¬ 
theses  in  Equation  (4.11).  This  definition  requires  the  introduction  of  some 
rather  general  fractal  arguments,  as  follows.  It  is  assumed  that  the  relation¬ 
ship  between  the  structural  character! sties  on  scales  m  and  s  can  be  repre¬ 
sented  by  a  single,  or  a  finite  number  of,  iterative  fractal  construction 
steps.  The  fractal  may  extend  to  scales  larger  than  s  or  smaller  than  m,  so  it 
represents  a  "covering"  or  "overlay"  onto  the  material's  structural  features  at 
those  scales.  Further,  a  different  fractal  may  be  employed  for  this  "overlay" 
for  each  state  of  the  system  (es,  cr$),  which  permits  the  ratio  of  characteris¬ 
tic  lengths  Ls/Lm  to  change  as  deformation  advances.  The  only  restriction  on 
the  fractal  is  that  the  dimensionalities  (D)  at  scales  m  and  s  be  equal,  so 
that  Equation  (4.11)  becomes 
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Next,  the  scale  factor  bm  is  defined  as 
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(4.13) 


where  the  characteristic  lengths  L$  and  Lm  represent  the  strained  condition  on 
each  scale,  consistent  with  the  above  description.  The  relation  =  N  is 
valid  for  any  fractal  description,  although  only  approximately  for  a  finite 
number  of  fractal  construction  iterations,  i.e.,  a  fractal  lattice.  Substitu¬ 
tion  of  this  relation  and  Equation  (4.13)  into  Equation  (4.12)  yields 
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where  Nm  is  defined  as  the  number  of  new  substructural  units  remaining  after 
the  fractal  iteration  (see  Section  2.2).  This  is  related  to  hardening  and 
damage  as  follows. 


The  creation  or  annihilation  of  hardening  or  damage  substructures  as 
deformation  advances  can  be  viewed  as  modifications  in  the  fractal  construction 
process.  Schematically,  this  is  shown  in  Figure  4.2.  On  the  scale  s,  the 
characteristic  length  L$  is  described  by  the  true  strain  e$ 

Ls  -  Lsq  e'5  (4.15) 


where  Lso  =  Ls(e$  =  o).  On  scale  m,  the  same  total  length  L$  is  made  up  of  Nm 
units  of  initial  length  Lso/bm,  which  are  each  strained  by 
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The  term  in  parentheses  is  lm,  analogous  to  L$  in  Equation  (4.15).  The  term  Nm 
can  be  represented  as  a  departure  from  the  homogeneous  condition,  which  is 
specified  by  Nm  =  bffl  =  3  in  Figure  4.2.  In  the  homogeneous  case,  the  line  on 
the  s  scale  remains  a  line  on  the  m  scale,  as  shown  by  the  dimensionality  Dm  = 
i.nNm/Jlnbm  =  *nbm/.fcnbm  =  1.  In  the  hardening  case,  extra  subunits  are  intro¬ 
duced,  so  that  Nm  >  bm  and  Dm  >  1.  The  three-sided  "box"  in  the  lower  left 
hand  side  of  Figure  4.2  represents  a  hardening  structure  such  as  a  dislocation 
tangle  or  a  cell  wall,  each  of  which  effectively  increases  the  magnitude  of  the 
energetic  path  for  the  passage  of  another  dislocation  from  point  A  to  B.  Con¬ 
versely,  the  open  section  at  scale  m  in  the  lower  right  hand  side  of  the  figure 
represents  a  void  or  other  similar  structure  which  decreases  the  resistance  for 
subsequent  deformation. 
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Defined  by  Simple  Fractal  Models,  with  bm  = 
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Nm  was  defined  simply  as 
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where  rm  >  1  for  domination  by  hardening  and  rm  <  1  for  domination  by  damage 
(softening) . 

4.1.3  Scale  Invariance  of  Energy  Storage  Rate 

Substitution  of  Equation  (4.17)  into  (4.14)  gives 
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The  derivative  on  the  right  hand  side  of  Equation  (4.18)  is  considered 
unknown,  a  condition  remedied  by  the  following  change  of  variables: 
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Insertion  of  this  equation  into  Equation  (4.18)  reveals  that  definition  of  the 


leaves  the  stored  energy  invariant  with  respect  to  the  strain  on  each  scale 
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This  is  the  key  relation  needed  for  development  of  the  stress  polarization 
equation.  Note  the  similarity  to  the  scale  invariant  stress  of  Section  3  under 
linear  elastic  conditions. 

4.1.4  The  Polarization  Equation 

Development  of  the  stress  polarization  equation  eventually  requires  that 
the  energy  storage  rate  terms  above  be  expressed  in  terms  of  effective  moduli 
so  that  the  path-dependent  DEM  formulations  in  [10]  can  be  employed.  This  is 
accomplished  as  follows.  For  any  sale,  the  energy  terms  (U)  can  be  expressed 
in  the  following  form,  where  subscripts  represent  principal  axes: 
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(4.22; 


where  \  is  Lame's  constant,  G  is  the  shear  modulus,  and  the  last  term  in  paren¬ 
theses  above  represents  the  sum  of  squares  of  shear  strains.  Because  Tre  =  0 
has  already  been  assumed,  di fferent i at i on  of  Equation  (4.22)  by  the  principal 
strains  t  gives 


dU.  dG.  7 

_L  =  2G.  e.  +  -a—  e. 
de.  i  i  d e.  l 

i  i 


(4.23) 


for  any  scale  i.  Substitution  of  this  equation  into  Equation  (4.21)  and 
employing  the  relation  2G-j  e-j  =  a-j  /  ( 1  +  v-j )  gives 
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where  v  is  Poisson  s  ratio,  which  is  dependent  on  the  dimensionality  as 
follows  [36] 
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However,  =  Dm  =  0,  so  vs  =  vm  =  v.  Rearranging  Equation  (4.24)  gives 


But  by  Equation  (4.9) 
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Definition  of  the  stress  polarization  at  scale  m  as  (see  Figure  4.1) 
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and  substitution  of  the  last  two  equations  into  Equation  (4.26)  gives  the 
stress  polarization  equation 


•  (1*v)  (dej;  4  -  °s  ^  '"<VEs>  •  TT71) 


where  the  modulus  term  represents  the  rate  of  mi crostructural  change  on  scale 
m,  and  the  last  bracket  effectively  contains  only  macroscopic  data  because  v 
can  be  approximated  by  the  usual  macroscopic  constant  value.  This  occurs 


because  changes  in  D  are  relatively  small  over  most  scales  and  strain  ranges  of 
interest,  an  observation  made  during  the  work  described  in  the  Appendix. 

It  should  be  noted  that  the  spatial  dependence  of  the  internal  stress 
fluctuations  does  not  appear  in  the  above  polarization  equation.  This  depend¬ 
ence  has  beer  absorbed  into  the  scale  factor  b^,  defined  in  Equation  (4.13). 
Also  the  "pol ari zati on"  itself  in  Equation  (4.29)  represents  the  departure  from 
a  "true  continuum,"  where  the  stresses  are  equal  on  all  scales  of  interest,  as 
in  a  fluid.  Of  primary  interest  is  the  envelope  defined  by  the  dashed  lines  in 
Figure  4.1,  rather  than  the  details  of  the  waveform  plotted  versus  the  X-axis 
in  that  figure.  The  scale-dependent  envelope  of  stress  fluctuations,  depicted 
by  the  heavy  dashed  curves,  is  actually  expected  to  exhibit  a  periodic  nature 
representati ve  of  nonlocal  behavior  confined  to  specific  scale  regimes  by  sub¬ 
structures  of  specific  characteristic  sizes.  The  solution  for  this  envelope 
shape  as  a  function  of  scale  will  yield  useful  information  concerning  substruc¬ 
ture  evolution,  and  is  the  subject  of  future  work. 

4.1.5  Path-Dependent  Effective  Moduli 

In  Equation  (4.29),  the  rate  of  change  of  the  effective  modulus  with 
respect  to  strain  is  given  by  the  path-dependent  differential  effective  medium 
(DEM)  theory  of  [10].  Path  dependence  means  that  the  effective  modulus  at  any 
state  (es,  <js)  depends  on  the  entire  sequence  of  combinations  of  hardening  and 
damage  substructures  that  occurred  to  reach  that  state.  The  DEM  approach 
describes  the  rate  of  change  of  the  effective  modulus  in  terms  of  the  amount 
and  rate  of  addition  of  these  two  types  of  inhomogeneities  to  a  parent  mate¬ 
rial.  These  two  inhomogeneities  are  described  in  terms  of  their  volume  frac¬ 
tions,  herein  denoted  4>h  and  4^,  where  the  subscript  m  has  been  dropped  for 
efficiency  of  notation.  The  volume  fraction  of  the  parent  material  (the  "back¬ 
bone")  is  denoted  4>0,  and  it  is  required  that 
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The  rate  of  change  of  the  bulk  (K)  and  shear  (G)  moduli  with  respect  to  an 
independent  variable,  denoted  here  by  (*),  are  given  by: 
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where  4>  =  4>h  +  4>d.  The  initial  conditions  are  K(0)  =  KQ,  G(0)  =  GQ.  The  P  and  Q 
above  depend  on  the  shapes  of  the  substructures,  and  were  selected  as  follows. 
Hardening  structures  were  assumed  to  be  of  plate  or  disk  form,  analogous  to 
grain  boundaries,  dislocation  tangles,  or  cell  walls.  The  moduli  of  the 
hardening  structures  (Kh,  Gh)  were  thus  assumed  to  be  those  representing  the 
state  of  maximum  hardness  in  the  macroscopic  material.  The  damage  structures 
were  assumed  to  be  of  a  spherical  void  form,  representing  a  ductile  material 
such  as  aluminum.  The  disc  shape  would  have  introduced  singularities  into  the 
formulation  when  the  damage  moduli  (Kd,  Gd)  were  assumed  to  be  zero,  as 
expected  physically.  The  P  and  Q  for  hardening  substructures  are  [10] 
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and  for  damage  substructures  are 
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Preliminary  results  concerning  the  verification  of  the  above  scale-  and 
rate-dependent  stress  polarization  formulation  (Equations  (4.29)  to  (4.35))  are 
shown  in  the  next  subsection. 


4.2  PRELIMINARY  VERIFICATION  RESULTS 

Limited  verification  of  the  above  scale-  and  rate-dependent  stress  polari¬ 
zation  equation  was  obtained  for  the  case  of  pure  aluminum  undergoing  creep  at 
3.0  MPa  and  593°K  [37].  Because  of  the  complexity  of  Equations  (4.29)  to 
(4.35),  some  simplifications  were  required  for  this  effort  to  remain  within 
project  resources.  Consequently,  verification  was  obtained  only  for  hardening 
behavior.  The  general  approach  was  to  introduce  approximations  for  some  of  the 
data  and  terms  needed  for  the  complete  solution  of  Equation  (4.29).  This  equa¬ 
tion  was  then  rearranged  so  that  the  characteristic  sizes  of  the  controlling 
substructures  could  be  predicted  as  a  function  of  the  macroscopic  strain,  and 
the  results  compared  to  data.  This  was  accomplished  as  follows. 

4.2.1  Specialization  for  Hardening  Behavior 

The  following  change  of  variable  in  the  modulus  rate  term  of  Equa¬ 
tion  (4.29)  was  found  to  eliminate  dependence  on  leaving  only  dependence 
on  <t>h: 


dG  dGm  dr 
m  _  m  m 

ds  dr  d  e 
m  mm 
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It  was  assumed  that  rm  was  adequately  approximated  by 


r 
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1  * 


(4.37) 


so  that  with  4q  =  =  0, 


(4.38) 


4.13 


where  the  value  of  G  =  f(<t>h(e$))  at  the  current  hardened  state  was  approximated 

by  the  larger  of  the  Hashi n-Shtri kman  bounds  for  two-phase  composites  [10], 

consistent  with  domination  by  hardening  (<t>d=0)  and  the  use  of  the  strain  energy 

in  Equation  (4.22)  [30,31].  Again  note  the  dropping  of  the  subscript  m  from  <J>, 
★ 

Gh,  and  terms  for  ease  of  notation. 

A  representati on  for  drm/defT1  was  found  by  combining  Equations  (4.15)  and 
(4.16)  to  give 


e  =  e  -  £nr  (4.39) 

m  s  m 

which  was  then  differentiated  with  respect  to  and  combined  with 
Equation  (4.20)  to  yield 


dr 

—  =  r  (r  b  -1 ) 
de  m  m  m 
m 


(4.40) 


This  equation  represents  the  case  of  multiaxial  deformation  because  the 
vertical  segments  in  the  lower  left  diagram  of  Figure  4.2  are  allowed  to 
strain.  A  uniaxial  representation  is  also  possible  if  Equation  (4.16)  is 
rewritten  so  that  vertical  segments  do  not  strain: 

L^freEm)b^(W(N”-b™)  <4-4i> 


which  when  combined  with  Equation  (4.15)  gives 

z 

e=*n(es  +  l-rm) 
m  m 


(4.42) 


Combination  of  this  equation  with  Equation  (4.20)  then  gives  the  unaxial 
representation 


Substitution  of  Equation  (4.40)  or  (4.43)  into  Equation  (4.29)  then  gives 


r J  r  b  -  l)  !  e  -*n  r  _p  =  R 
m  m  m  ■  s  m 


(4.44 


for  the  multiaxial  case,  or 


e  5fr  b  -l)  +  (r  -llj  [*nfe  s+  1  -  r  ) J ^  =  R 
•mm  m  m 


(4.45 


for  the  uniaxial  case,  where 


l«s(»*”(VS)  -iw)  +T^J  (‘-*h) 


G  +  G. 

(G  .  G)  - 

Gh  +  Gh 


(4.46 


4.2.2  Calculations 

Equations  (4.44)  to  (4.46)  were  used  to  solve  for  the  strain-dependent 
evolution  of  the  controlling  microstructure  sizes,  as  follows.  Strain  rate 
versus  strain  data  for  the  creep  of  pure  aluminum  at  a  shear  stress  of  3.0  MPa 
and  593°K  were  taken  from  Figure  2a  of  [37],  There  were  no  strain  data  for 
<  0.2,  so  the  initial  value  for  computation  was  assumed  to  be  0.002,  the 
point  at  which  plastic  deformation  is  generally  considered  to  be  fully 
developed.  The  strain  history  was  thus  es  =  0.002,  0.2,  0.3,  0.4,  0.5,  G.6, 
and  the  corresponding  strain  rate  history  was  e  =  5  x  10  ,  4.5  x  10  ,  1  x 


10  ,  7  x  10“  ,  6  x  10"  ,  5.5  x  10“b  (sec-1).  Also,  eq  and  '?  were  assumed  to 
be  10“2/sec  and  1/3.1,  respectively.  The  material  constant  C  was  again  taken 
to  be  unity.  The  initial  value  of  the  shear  modulus  (G0)  was  taken  as 


mmmfssmsm 


13.45  GPa,  with  Poisson's  ratio  v  =  0.347.  The  remaining  unknowns  were  then 
Gh,  $h,  and  op  as  functions  of  es.  The  stress  polarization  history  was  assumed 


to  be  op  =  0.1,  0.1,  0.2,  0.3,  0.4,  0.5.  There  were  no  data  concerning  the 


hardening  fraction  evolution  with  strain,  although  inspection  of  the  micro¬ 
graphs  in  [37]  resulted  in  the  estimate  =  0.3  at  =  0.6.  The  assumed 
hardening  fraction  history  was  =  0.025,  0.15,  0.2,  0.25,  0.275,  0.3.  was 
assumed  to  be  slightly  larger  than  GQj  and  is  discussed  in  a  later  paragraph. 

The  above  data  were  input  to  a  small  computer  code  that  solved  numerically 
for  the  scale  factor  bm  =  f(es)  at  which  the  substructure  controlled  the 
macroscopic  strain  rate.  This  bm  was  then  used  to  calculate  the  evolution  of 
the  controlling  substructure  size  Lm  from  Equation  (4.13)  with  the  reference 
length  Lso  =  1  cm: 


L„  * 


(4.47) 


Results  for  both  the  uniaxial  and  multiaxial  cases  are  shown  in  Fig¬ 
ure  4.3,  and  agree  reasonably  well  with  available  data,  despite  the  large 
number  of  assumptions  required  to  complete  these  preliminary  calculations.  At 
small  strains,  the  controlling  substructure  size  is  close  to  the  dislocation 
spacing  for  the  estimated  initial  dislocation  density  of  106/cm^,  a  reasonable 
value  for  pure  aluminum  in  the  annealed  state.  The  controlling  substructure 
size  (dislocation  spacing)  decreases  rapidly  as  plastic  deformation  advances 
and  the  dislocation  density  increases  drastically.  There  are  no  met a  11 ographi c 
data  for  es  <  0.2,  and  the  shaded  area  in  Figure  4.3  indicates  the  anticipated 
range  possible  for  randomly  distributed  peak  dislocation  densities  between  10^ 
and  1010/cm^.  At  es  =  0.05  to  0.1,  dislocation  tangles  would  be  expected  to 
begin  forming  the  walls  of  cells  or  subgrains.  These  walls  would  decrease  in 
width  as  deformation  advances,  thus  increasing  the  available  path  length  for 
the  motion  of  mobile  dislocations  within  the  cells  or  subgrains.  This  is  indi¬ 
cated  by  the  general  upward  trend  in  Figure  4.3,  which  begins  to  stabilize  at 
larger  strains,  and  results  in  a  good  prediction  for  the  subgrain  sizes  found 
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FIGURE  4.3.  Predicted  Evolution  of  Controlling  Substructure 
for  Creep  of  Aluminum  at  3  MPa,  593°K 


4.2.3  Relation  to  Physics  and  Metallur 


The  above  results  were  obtained  by  assuming  particular  values  for  the 
ratio  Gh/GQ,  as  shown  i n  Figure  4.3.  These  values  were  found  to  be  reasonable 
by  the  following  arguments.  GQ  is  the  shear  modulus  of  pure  aluminum  in  the 
initial  annealed  state,  while  G^  is  the  modulus  of  the  hardened  phase  (i.e.,  a 
cell  wall)  where  dislocation  densities  are  large  (about  1.5  x  lO^/cm2  [38,39]) 
and  atomic  bonds  are  elastically  strained  in  an  average  sense.  The  force- 
displacement  curve  of  atomic  theory  exhibits  an  increase  in  slope  as  the  atoms 
are  slightly  separated  [40],  This  increase  in  interatomic  force  with  small 
strains  is  manifested  as  an  apparent  increase  in  the  effective  modulus,  which 
is  defined  as  the  ratio  of  stress  (force  per  unit  area)  over  strain.  The  force 
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increases  in  Figure  4.3  are  6  x  10  to  1.25  x  10"'>  percent,  and  may  not  be 
easily  measurable.  However,  analysis  of  the  corresponding  increase  in  strains 
is  possible,  as  follows. 

A  “universal  equation  of  state  for  metals"  is  available  from  [40],  where 
the  interatomic  potential  E  is  expressed  as 


E  =  E  E*(a*) 
o  ' 


(4.48) 


where  E0  is  the  equilibrium  binding  energy  and 


a*  =  (r-rg  )/ 1 


(4.49) 


with  r  as  the  Wigner-Seitz  radius,  rg  =  1.58  x  10-8  cm  is  the  equilibrium  (zero 
strain)  value  of  r  for  aluminum,  and  l  =  0.336  x  1U'8  cm  is  the  normalization 
factor  for  aluminum.  The  normalized  potential  is 


E*  =  e  (1  +  a*  +  0.05  a*J) 


The  tangent  modulus  is  proportional  to  Ao/Ae, 


r  do  1  dF 
b  =  dl  =  T  37 

where  the  constant  area  A  contains  the  Poisson  factor  and  F  is  force, 
strain  is  defined  as 


e  =  a  *1/1 


where  £0  =  4.04  x  10-8  cm  is  the  lattice  constant.  Thus, 


(4.50) 


(4.51) 


(4.52) 


AX/  cla^ 


(4.53) 
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(4.54) 


Gh  (d2E*/da*2)h 
Go  (d2E*/da*2)Q 

so  that  E0  and  the  constants  in  Equation  (4.53)  cancel.  Taking  the  required 
derivatives  and  assuming  that  a^  >  a*  =  o,  we  have 

ji  =  e"a*  (1  +  2.7a*  -  a*2  -  0.65a*3  -  0.05a*5)  (4.55) 

o 

Setting  Equation  (4.55)  equal  to  1.000006  or  1.0000125  yields  a*  =  3.5  x  10-6 
or  7.4  x  10"6  for  the  uniaxial  and  multiaxial  cases  in  Figure  4.3,  respec¬ 
tively.  Equation  (4.52)  then  gives  the  average  elastic  strains  in  the  hardened 

phase;  eh  =  2.9  x  10-7  (uniaxial),  =  6.2  x  10-7  (multiaxial).  These  values 

may  be  compared  to  the  average  strains  associated  with  dislocations,  as 
fol lows. 

The  volume  fraction  of  dislocations  in  the  hardened  phase  may  be  estimated 
as 


Vf  =  *Uo)2  P  =  7.7  x  10"6  (4.56) 

Assuming  that  the  material  between  the  dislocations  (whose  spacing  is  p"^7^  = 
643  Aq)  has  zero  elastic  strain,  and  considering  the  dislocations  as  "fibers" 
so  that  the  rule  of  mixtures  applies,  the  strain  associated  with  the  disloca¬ 
tions  is  e^/Vf  =  0.037  (uniaxial)  or  0.081  (multiaxial).  The  shear  strain 
associated  with  a  screw  dislocation  is  b/2uJl0  =  0.113  (b  =  2.86  x  10'^  cm  is 
Burger's  vector).  The  strain  associated  with  a  path  around  an  edge  dislocation 
is  b/8£Q  =  0.088.  Therefore,  the  strains  computed  from  the  analysis  of  the 
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modulus  changes  (4  to  8%)  compare  reasonably  well  with  those  from  a  metallurgi¬ 
cal  analysis  (9  to  11%).  The  factor  of  two  discrepancy  is  easily  accounted  for 
by  the  uncertainties  in  the  estimated  dislocation  density,  and  the  approxima¬ 
tion  of  a  constant  value  for  over  all  scales.  This  lends  support  to  the 
verification  of  the  stress  polarization  equation. 
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4.3  DISCUSSION 

The  scale-  and  rate-dependent  stress  polarization  ( ap )  equation  was 
developed  from  Bodner’s  Equation  using  fractal  scaling  relations  in  a  manner 
similar  to  the  mi crostructural  form  developed  in  Section  3.  The  latter  form 
(Equation  3.12)  assumes  that  local  behavior  dominates  with  respect  to  scale, 
and  requires  mi crostructural  models  to  predict  the  macroscopic  strain  rate. 
Conversely,  the  ap  formulation  (Equation  4.29)  quantifies  the  nonlocal  behavior 
as  a  function  of  scale.  It  requires  no  mi crostructural  models,  but  needs 
certain  state  variable  histories  to  predict  the  evolution  of  the  internal 
stress  field  fluctuations  as  a  function  of  the  macroscopic  stress  and  strain 
rate. 

The  preliminary  verification  results  for  aluminum  (Section  4.2)  showed 
that  the  ap  equation  could  be  used  to  predict  the  evolution  of  the  micro- 
structural  features  that  control  viscoplastic  deformation,  given  the  hardening 
and  damage  state  variable  histories.  Although  such  data  are  sometimes  avail¬ 
able  for  larger  strains,  they  frequently  do  not  exist  at  smaller  strains  (where 
viscoplastic  rate  effects  are  most  important)  because  of  experimental  diffi¬ 
culties.  A  method  has  been  developed  to  alleviate  this  deficiency  by  eliminat¬ 
ing  ^  and  4>dm  from  the  ap  equation,  leaving  the  stress  polarization  as  a 
function  of  the  chosen  scale  factor  (bm)  and  the  macroscopic  stress  and  strain 
rate.  The  resulting  equation  can  then  be  solved  to  find  the  spectrum  of  scales 
over  which  nonlocal  behavior  dominates,  thus  providing  guidance  for  metal- 
lographic  analyses  and  possibly  for  metallurgical  design. 

Note  that  the  ap  equation  contains  the  function  rm  =  f ( 4>hm , 4>dm ) »  which  has 
two  important  features.  First,  concerning  the  verification  results  in 
Figure  4.3,  it  should  be  noted  that  [37]  contains  mention  of  a  small  fraction 
of  the  subgrains  having  cells  with  average  diameters  an  order  of  magnitude 


smaller  than  the  subgrains.  References  [38,  39]  also  found  substructure  domi¬ 
nated  by  cells  about  8  microns  in  diameter.  It  is  suspected  that  more  accurate 
and  complete  representations  of  rm,  and  thus  would  be  more  complicated 
mathematically,  resulting  in  multiple  solutions  for  bm  (branchings)  that  would 
predict  both  cells  and  subgrains.  Second,  rm  measures  hardening  and  damage  in 
terms  of  departures  from  unity.  Expansion  to  multiaxial  forms  can  provide  a 
combined  hardening/damage  tensor  that  can  be  inserted  into  continuum  formu¬ 
lations,  thus  completing  the  second  objective  in  the  first  paragraph  of  the 
Summary . 
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5.0  CONCLUSIONS  AND  RECOMMENDATIONS 


The  following  results  were  obtained  during  the  course  of  this  work: 

•  A  new  form  of  Bodner's  viscoplastic  constitutive  equation  was  developed. 
This  form  (Equation  1.1  or  3.12)  admits  mi crostructural  models  in  an 
organized  and  physically  meaningful  manner.  It  was  verified  by  comparison 
to  high  temperature  creep  data  for  a  steel  alloy,  and  is  expected  to 
exhibit  improved  hardening  behavior  in  computer  simulations  of  structural 
performance .  The  origin  of  the  empi ri cal ly-determi ned  stress  exponent  was 
also  identified  in  terms  of  mi crostructural  phenomena  occurring  at  various 
characteri sti c  length  scales  within  a  deforming  material. 

A  stress  polarization  equation  for  viscoplasticity  was  also  developed  from 
Bodner's  Equation.  This  new  formulation  (Equation  1.2  or  4.29)  describes 
the  rate-dependent  evolution  of  the  internal  stress  field  fluctuations  in 
the  material  as  a  function  of  macroscopic  stress  and  strain  rate,  two 
hardening  and  damage  state  variables,  and  the  characteristic  sizes  of  the 
rate-controlling  microstructures.  Preliminary  verification  was  obtained 
for  the  creep  of  aluminum,  and  showed  that  the  stress  polarization  equa¬ 
tion  is  potentially  useful  for  the  prediction  of  mi crostructural  evolution 
in  terms  of  macroscopic  stress-strain  rate  data.  This  equation  also  con¬ 
tains  the  basis  for  the  development  of  a  new  material  state  tensor  that 
simultaneously  describes  hardening  and  damage  phenomena. 

The  following  recommendations  are  also  made: 

•  The  verification  results  obtained  herein  are  preliminary  in  nature,  and 
further  testing  of  the  above  equations  is  in  order. 

•  Development  of  the  stress  polarization  equation  should  be  continued  to 
(1)  provide  analytical  tools  capable  of  predicting  mi crostructural 
evolution  in  terms  of  macroscopic  response,  thus  giving  guidance  to 
metal  1 ographic  analyses  and  reducing  the  costs  of  metallurgical  design, 
and  (2)  complete  the  development  of  the  combined  hardening  and  damage 
state  tensor. 
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THE  CONSERVATION  OF  MASS  APPROACH 


This  approach  was  an  attempt  to  solve  for  the  dimensionalities  in  Equa¬ 
tion  (4.1),  as  follows.  The  conservation  of  mass  equation  requires  that  the 
material  derivative  of  the  mass  is  zero.  Expressing  the  mass  m  as  the  density 
times  the  fractal  volume,  this  is 


dm 

dt 


3  f 
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3  r  ,  D  s 
Vi  3x^  ^  ^ 


(A.l) 


where  the  small  d  represents  the  material  derivative  and  the  summation  conven¬ 
tion  applies  in  the  last  term. 

Performing  the  derivatives,  this  equation  becomes 
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The  usual  conservation  of  mass  equation  gives 
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which  can  be  substituted  into  Equation  (A. 2)  to  give  a  differential  equation 
for  the  fractal  volume  (L^): 


5L°  + 
W~  + 


3L°  .  D  dvi 

vi  3xT  "  L 


(A.4) 


I 


At 


,V.vL>: 


i.Xv'Cv.v'I-/ 


where  the  L  above  is  a  suitable  average  of  the  characteristic  lengths  of  the 
three  principal  directions 
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where  L0  is  the  initial  length,  and  v  is  Poisson's  ratio.  Equations  (A. 4)  and 
(A. 5)  are  assumed  to  be  valid  for  any  scale  within  the  material.  The  x^  and  v^ 
in  Equation  (A. 4)  are  defined  by  Equations  (A. 5)  and  their  derivatives. 

The  choice  of  the  average  length  L  determines  the  final  form  of  the  dif¬ 
ferential  equation  for  D.  A  number  of  attempts  were  made  to  find  a  suitable 
representati on  of  L,  with  inconsistent  results.  Examples  are  choices  based  on 
the  strain  invariants  (i.e.,  I^)  and  the  strain  deviator  invariants  (i.e., 

IV): 
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The  final  differential  equation  for  D  on  any  scale  corresponding  to,  i.e., 
Equation  (A. 7)  was 
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where  Poisson's  ratio  v  depends  on  the  dimensionality  D  [36] 


6  -  0 
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(A. 11) 


The  differential  equation  corresponding  to  Equation  (A. 6)  was  similar,  but  with 
different  versions  of  Equations  (A. 10). 

Given  a  characteristic  length  LQ  corresponding  to  a  chosen  scale,  dif¬ 
ferential  equations  such  as  Equation  (A. 8)  were  solved  numerically  and  the 
results  represented  in  the  form  of  plots  of  macroscopic  strain  rate  versus 
strain.  This  was  done  by  approximating  the  constitutive  equation  by  the  form 
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(A. 12) 


where  aQ  is  a  lattice  constant,  bQ  is  a  Burger's  vector,  ¥  is  the  Schmid 
factor,  EQ  is  atomic  binding  energy,  and  is  the  mobile  dislocation  density. 
The  term  f(0a)  represents  damage  on  the  atomic  scale,  and  is  developed  from  a 
rearrangement  of  Equation  (A. 8)  so  that 
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Note  that  l_s°s  is  omitted  from  Equation  (A. 12)  because  it  is  assumed  relatively 
constant,  consistent  with  the  observation  of  volume  conservation  on  the  macro¬ 
scopic  scale  [41],  Setting  the  last  term  in  Equation  (A. 12)  equal  to  a  con¬ 
stant,  the  basic  strain  rate  versus  strain  behavior  can  be  observed  as  shown  in 
Figure  A.l.  The  ordinate  in  this  figure  is  the  macroscopic  strain  rate  normal¬ 
ized  to  the  value  at  the  assumed  starting  strain  (10"^).  The  abscissa  is  the 
strain  on  the  scales  noted  in  the  figure. 

The  two  solutions  in  Figure  A.l  show  quite  different  results  depending  on 
the  representation  of  L  chosen.  Another  representati on  with  1^=0  in  Equa¬ 
tion  (A. 6)  [41]  gave  results  similar  to  the  "strain  invariant"  curves  shown  in 
the  figure.  Note  that  both  solutions  show  macroscopic  effects  at  much  smaller 
strains  than  for  the  microscopic  scale.  This  may  be  somewhat  consistent  with 
the  following  observation.  Although  a  small  permanent  deformation  of  the 
macroscopic  sample  may  be  measurable,  a  random  microscopic  sample  will  fre¬ 
quently  show  no  obvious  evidence  of  such  deformation,  at  least  for  the  rela¬ 
tively  small  macroscopic  strains  over  which  viscoplasticity  is  important. 
However,  the  monotonic  behavior  of  either  solution  in  Figure  A.l  was  taken  as 
an  indication  that  hardening  and  damage  were  actually  not  simultaneously  repre¬ 
sented,  as  originally  intended.  This  is  also  obvious  from  Equation  (A. 12), 
where  the  damage  term  may  be  represented  by  the  first  bracket  and  the  hardening 
behavior  by  pm  in  the  second  bracket.  In  essence,  the  continuum  and  metal¬ 
lurgical  contributions  had  remained  separated.  The  consequent  failure  of  the 
above  conservation  of  mass  approach  thus  indicated  that  another  direction 
should  be  taken,  and  the  stress  polarization  approach  was  subsequently  adopted. 
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